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Abstract 

We perform a numerical investigation of the field evolution in the 
smooth hybrid inflationary model. We find that for almost all the 
examined initial values we do get an adequate amount of inflation. 
Our results show that the model is "natural" and satisfactory. 
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The hybrid inflationary scenario^ 1 -* proposed by Linde in the context of 
non supersymmetric (SUSY) theories is a realization of chaotic inflation 
based on a coupled system of two scalar fields one of which may not be 
a gauge singlet. The great advantage of this scenario is that it produces the 
observed temperature fluctuations in the cosmic background radiation (CBR) 
with natural values of the coupling constants. However, inflation terminates 
abruptly and is followed by a "waterfall" regime during which topological 
defects can be easily produced. Recently, two of us have proposed^ a vari- 
ant of Linde's potential which can be derived in a wide class of SUSY grand 
unified theories (GUTs) based on semi-simple gauge groups by utilizing the 
first non-renormalizable contribution to the superpotential. Although one 
gets only a slight variation of Linde's potential, the cosmological scenario 
obtained is drastically different. Already since the beginning of inflation, 
the system follows a particular valley of minima which leads to a particular 
point of the vacuum manifold. Thus, this inflationary scenario does not lead 
to production of topological defects. Also the termination of inflation is not 
as abrupt as in the hybrid case. It is quite smooth and resembles more the 
cases of new or chaotic inflation. The main advantage of this smooth hybrid 
inflationary scenario is that the measured value of the temperature fluctua- 
tions of CBR can be reproduced with natural values of the parameters and 
with a GUT scale, Mx, consistent with the unification of the minimal super- 
symmetric standard model (MSSM) gauge couplings. It is also remarkable 
that the scale controlling the non-renormalizable terms in the superpotential 
turns out to be of order 10 18 GeV. The spectral index of density fluctuations 
is close to unity. 
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For an inflationary scenario to be considered fully successful, one has to 
show that it is obtainable for a wide class of "natural" initial values of the 
fields and their time derivatives. In ref.(2), two sets of initial conditions were 
studied semianalytically and it was argued that they lead to smooth hybrid 
inflation. These sets, however, cannot be considered completely "natural" 
since they require some or even a considerable discrepancy between the ini- 
tial values of the fields. The purpose of this paper is to identify a wide class 
of "natural" initial conditions for smooth hybrid inflation, i.e., comparable 
initial values of the fields for which the system falls at the bottom of a par- 
ticular valley of minima. Its subsequent evolution along this valley, then, 
produces smooth hybrid inflation. To this end, we solve numerically the evo- 
lution equations of the system for a wide class of initial conditions. The result 
is striking and unexpected. We find that, for almost all the examined initial 
conditions (except a narrow transition region), we do get smooth hybrid in- 
flation with an adequate number of e-foldings. This result together with the 
other advantages of this inflationary scenario makes it fully satisfying and 
"natural". Our analysis also includes cases with all initial field values being 
much smaller than the Planck scale. In these cases, our results are expected 
to be less affected by replacing global by local supersymmetry. 

The smooth hybrid inflationary scenario can be realized in the context 
of a SUSY GUT based on a gauge group G of rank > 5. We assume that 
G breaks spontaneously directly to the standard model (SM) group Gs at 
a scale M x ~ 10 16 GeV and that below M x the only SM non-singlet states 
of the theory are the usual MSSM states. This guarantees the successful 
MSSM predictions for sin 2 9 w and a s . The theory could also possess some 
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global symmetries. The symmetry breaking of G to Gs is obtained through 
a superpotential which includes the terms 

Here ,0 is a conjugate pair of left-handed SM singlet superfields which 
belong to non-trivial representations of the gauge group G and reduce its 
rank by their vacuum expectation values (vevs), s is a gauge singlet left- 
handed superfield, \i is a superheavy mass scale related to M x , whereas 
M is a mass scale of the order of the "compactification" scale M c ~ 10 18 
GeV which controls the non-renormalizable terms in the superpotential of 
the theory. The superpotential terms in eq.(l) are the dominant couplings 
involving the superfields s ,0 ,0 consistent with a continuous R-symmetry 
under which W — > e te W,s e i9 s, 00 — ■> 00 and a discrete symmetry under 
which 00 changes sign. 

The potential obtained from W in eq.(l), in the supersymmetric limit, is 

V = |p» - <f!| 2 + 4|.|»1^(|#|' + + D - terms, (2) 

where the scalar components of the superfields are denoted by the same sym- 
bols as the corresponding superfields. Vanishing of the D-terms is achieved 
along the D-flat directions where |0| = |0|. The supersymmetric vacuum 

< s >= 0, < >< >= ±[xM, | < > | = | < > | (3) 

lies on the particular D-flat direction 0* = ±0. Restricting ourselves to this 
direction and performing appropriate gauge, discrete and R-transformations 
we can bring the complex s,0,0 fields on the real axis, i.e., s = = = 
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\Xi where a and x are re& l scalar fields. The potential in eq.(2) then takes 
the form 

and the supersymmetric minima correspond to 

| < X > I = 2(/iM) 1/2 , < a >= . (5) 

The mass acquired by the gauge bosons is Mx = g(/iM) 1//2 , where g is the 
GUT gauge coupling. For any fixed value of a, the potential in eq.(4), as a 
function of x 2 , has a local maximum at x 2 — an d an absolute minimum 
lying at 



X 2 --^,/ora 2 »/i¥, (6) 



2 _ 4^ 
" ~' 3 a 2 

and x 2 — 4/iM, for <x 2 << /jM. The value of the potential along the maxima 
at x 2 = is constant, V max (x 2 = 0) = 

Assume for the moment that at some region of the universe the scalar 
fields x an d cr, starting from appropriate initial conditions to be discussed 
below, evolve in such a way so they become almost uniform with values at 
the bottom of the valley of minima in eq.(6). One can then show^ 2 ) that 
these fields, at subsequent times, move towards the supersymmetric minima 
in eq.(5) following this valley and the system inflates till a reaches the value 

CT — ( v^ )1/S( " M) ' /2 1 ( " M) ' /2 ' (7) 

where Mp = 1.2 x 10 19 GeV is the Planck mass. The number of e-foldings 
from the moment at which the a field has the value a till the end of inflation 
is given by 
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This implies that the value of the a field when the present horizon crossed 
outside the inflationary horizon was 

^M^o, 0) 

where Njj is the number of e-foldings of the present horizon size during in- 
flation. After the end of inflation, the a and \ fields enter smoothly into an 
oscillatory phase about the global supersymmetric minimum of the potential 
in eq.(5) with frequency m a = m x = 2v / 2(/i/M) 1 / 2 /i. These fields should 
eventually decay into lighter particles and "reheat" the universe. Taking 
N H = 60, M x = 2 x 10 16 GeV, g = 0.7 (consistent with the MSSM unifica- 
tion), the microwave background quadrupole anisotropy (AT/T) ~5x 10~ 6 
from the Cosmic Background Explorer (COBE) implies that M ~ 9.4 x 10 17 
GeV and fi ~ 8.7 x 10 14 GeV. With these numbers, which we will use through- 
out this paper, one estimates the value of a at the end of inflation (see eq.(7)) 
a Q ~ 1.1 x 10 17 GeV ~ 5AMx and its value when the present horizon size 
crossed outside the inflationary horizon (see eq.(9)) a H — 2.54a ~ 2.7 x 10 17 
GeV. An important advantage of the above smooth hybrid inflationary sce- 
nario is that inflation takes place at relatively low values of the field a and, 
therefore, there is hope that it survives even in the context of supergravity 
theories although it may acquire drastic modifications. 

We will now try to specify the initial conditions for the a and x fields 
which lead to the above described inflationary scenario. In other words, we 
will try to identify the initial conditions for which the system falls at the 
bottom of the valley of minima in eq.(6) with a value of a > <jh so that its 
subsequent evolution along the valley produces an adequate amount of infla- 
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tion. We assume that, after "compactification" at some initial cosmic time, 
a region emerges in the universe where the scalar fields a and x happen to be 
almost uniform with negligible kinetic energies. (The initial values of a and 
X can always be transformed by appropriate gauge and R-transformations to 
become positive.) The evolution of the system in this region is governed by 
the following equations of motion 

* + 3ff *-2^ 2 -W> + W = °' (10) 

where overdots denote derivatives with respect to cosmic time and H is the 
Hubble parameter 

H = )^ Mp Y /2 = (f f /2 M P \\x 2 + \* + V( X , <x)) 1/2 (12) 

(g is the energy density). 

The case where initially a ^> Mp ^> x was examined in ref.(2).Under 
these circumstances, the last term in eq. (4) is initially the dominant contribu- 
tion to the potential energy density of the system (assuming x ~ (/U 2 M 2 /<r) 1 / 3 ). 
Also, eq.(10) reduces to 

* +3 ^ + w^°- (13) 

Let us for the moment assume that a remains almost constant. The frequency 
of oscillations of the x field is then much greater that H and x initially 
performs damped oscillations over the maximum at x — 0- The continuity 
equation 

q = -3H(q + P ), (14) 



with p being the pressure averaged over one oscillation of x becomes 



q=-ZH 16 , (15) 

where 7 = 3/2 for a x 6 potential 1 ^. This equation together with the fact 
that q is proportional to H 2 gives 

4 

H~—. (16) 
9t v 1 

Comparing this result with eq.(12) we can obtain the amplitude of the oscil- 
lating x held 

Eq.(ll) averaged over one oscillation of x then gives 

The solutions of this equation are of the form a = t a , where a satisfies the 
quadratic equation 

, 1 1 ,M P ^ , x 

a+ r + 2fi^^°- (19) 

For a 3> Mp, the solutions are 

1 1 ,M Pn9 , 1 ,M Pn9 , . 

a~ — + — — anda- — • 20 

3 9tt a 1 9vr a 1 v 1 

This means that a quickly approaches an extremely slowly decreasing func- 
tion of time and, thus, our starting assumption that it remains approx- 
imately constant is justified. When the amplitude of the x field drops 
to about (/i 2 M 2 /(r) 1 / 3 , the /i 4 term dominates the potential in eq.(4) and 
the Hubble parameter becomes approximately constant and equal to H = 
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(87r/3) 1//2 /i 2 /Mp and remains so thereafter till the end of inflation. The sub- 
sequent evolution of the system has been studied in detail in ref.(2). The 
overall conclusion is that, in a time interval 



the x held falls into the valley of minima in eq.(6) and relaxes at the bottom of 
this valley whereas the a field still remains unchanged and much greater than 
Mp. After that, the system follows the valley of minima towards the super- 
symmetric vacuum and, therefore, the smooth hybrid inflationary scenario is 
realized for initial values of the fields satisfying the inequality a 3> M P 3> X- 
However,these initial conditions cannot be considered totally satisfying be- 
cause, assuming that the initial energy density is well below M P , we see that 
there must be some discrepancy between the initial values of the fields. Also 
the inclusion of supergravity is expected to invalidate the above discussion of 
initial conditions which involve values of the field a > Mp. To minimize the 
influence from supergravity one could start with field values much smaller 
than M P and, as pointed out in ref.(2),still obtain adequate inflation. How- 
ever, this required an initial energy density much smaller than Mp and an 
initial value of x "unnaturally" smaller than the initial value of a. 

For smooth hybrid inflation to be considered as a fully successful infla- 
tionary scenario, one must show that it is obtained for a wide class of initial 
conditions which are more "natural" than the ones just discussed. This can 
be done only numerically. To this end, we put x = x/Mp and a = a /Mp in 
eqs. (10) and (11) which become 



At ~ Qk(^-) 2 H- 1 




(21) 



X " + SH X ' - -^-(16/z 2 M 2 - x 4 ) + 
32M 4 



3x 5 a 2 



= 0, 



(22) 



8M 4 
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where 



with 



a" + 3H&' + = , (23) 
8M 4 



H = J- = (y ) 1/2 [i(x') 2 + + VU -)] 1/2 (24) 



*) - %^ = (A 2 " ^) 2 + ■ (25) 

y ' Mj, VP 16M 2; 16M 4 V ; 

Here primes denote derivatives with respect to the dimensionless time vari- 
able r = M P t and 

M 1 , /i 1 

M = ~ a = ~ (26) 

M P 12.83' P M P 13737' v ; 

We have conducted numerical integration of eqs.(22) and (23) for an extensive 
set of initial values of the fields in the ranges 0.01 < a < 1.2 and 0.01 < 
X < 0.5 and with vanishing initial velocities. The integration of these two 
coupled equations was performed by implementing a variant of the Bulrish- 
Stoer^ 4 ) variable step method in a Fortran program run mainly on a number 
of workstations. As a general rule, the initial step for the dimensionless time 
variable r was chosen to be 10 while the sought accuracy was put to 10~ 12 . 
This choice was found to ensure reasonable stability in the cases where a 
settles down to a constant non-zero value relatively early. In the opposite 
cases as well as near some transition regions, the initial step was decreased 
to 5 or less, while the sought accuracy was increased by two or three orders 
of magnitude depending on the case. 

The results of our search are summarized in figures 1 and 2 . Each point 
shown on the a-% plane corresponds to a given set of initial conditions and 
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depicts a definite evolution pattern for the a- x system according to the sym- 
bol attributes being used. We have used filled circles to specify the evolution 
pattern where both fields oscillate and fall rapidly to the supersymmetric 
minima in eq.(5) without producing any appreciable amount of inflation. 
Open triangles correspond to the case where a starts-off at relatively large 
values (cr > x) an d decreases slowly tending asymptotically to a constant 
value a >Bh ~ 0.0225 . The field x oscillates and relaxes at the bottom of 
the valley in eq.(6). The system then evolves through the valley of minima 
in eq.(6) giving an adequate amount of inflation. Finally, open circles cor- 
respond to the pattern where both fields start oscillating at the beginning. 
Then, a settles down at large values and the system subsequently follows the 
valley of minima in eq.(6) as in the preceding case. 

The evolution pattern represented by open triangles includes the pat- 
tern found in the limiting case a 1 ^> x analyzed earlier by means of 
semianalytic arguments. In general, however, a does not remain frozen as 
in the limiting case, but its variation over a large period of time is small. 
The open circles area being the least expected deserves further attention. 
Here, although a starts at moderate or small values, it appears to increase 
in amplitude absorbing energy from the fast oscillating field x an d creates 
conditions that eventually lead to evolution of the open triangles type. It 
is a beautiful example of large energy transfer between two strongly-coupled 
non-linear oscillators. 

The points depicted in figs.l and 2 follow a remarkably regular pattern, 
although there is some intermingling. This intermingling persisted even when 
the sought accuracy was increased to a maximum and the step was lowered 
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to a minimum, thus, we have to assume that it really exists. It would, of 
course, be very desirable to get three distinct regions separated by critical lines 
but we have no reason to a priori exclude the possibility of intermingling. 

The problem to be addressed next concerns the "naturality" of initial con- 
ditions. A closer look at fig. 1 reveals that the open triangles area,although 
leading to successful inflation, cannot be considered as being completely "nat- 
ural" since it requires relatively large differences between the initial values of 
the fields. The open circles region appears to be significantly better since the 
initial values can be of the same order of magnitude. Considering that the 
initial kinetic energies for both fields are taken to be zero, the initial energy 
density equals the potential energy density given in eq.(25). It turns out 
that, when x an d <? are of the same order of magnitude and approximately 
equal to a few tenths, we get acceptable initial energy densities. Thus, part 
of the open circles area corresponds to "natural" initial conditions which 
lead to successful inflation. Inclusion of supergravity will certainly invalidate 
the preceding analysis, except for an area where all initial field values are 
much smaller that the Planck scale. This area is shown in fig. 2 and a closer 
inspection shows that all but one of the points fulfil all the conditions for ad- 
equate inflation. In summary, taking into account all the previously stated 
results, we feel confident to conclude that the smooth hybrid inflationary 
model appears to be "natural" and satisfactory. 
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Figure Captions 



Fig. 1: Evolution patterns for the <7— X system. Filled circles represent points that 
do not lead to inflation. Open triangles give adequate inflation having only the X ne ld 
oscillating. Open circles give adequate inflation with both fields oscillating initially. 

Fig. 2: Same as in Fig.l. The initial values are now restricted to lie near the beginning 
of the axes. 
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